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Abstract. It is shown that one can count fc-edge paths in an n-vertex 
graph and m-set fc-packings on an n-element universe, respectively, in 
time (^."2) and {^f.^, up to a factor polynomial in n, k, and m; in 

polynomial space, the bounds hold if multiplied by 
spectively. These are implications of a more general result: given two set 
families on an n-element universe, one can count the disjoint pairs of sets 
in the Cartesian product of the two families with 0(nt) basic operations, 
where £ is the number of members in the two families and their subsets. 



1 Introduction 

Some combinatorial structures can be viewed as two halves that meet in the 
middle. For example, a fc-edge path is a combination of two fc/2-edge paths. 
Bidirectional search [10,25] finds such structures by searching the two halves 
simultaneously until the two search frontiers meet. In instantiations of this idea, 
it is crucial to efficiently join the two frontiers to obtain a valid or optimal 
solution. For instance, the meet-in-the-middle algorithm for the Subset Sum 
problem, by Horowitz and Sahni [15], implements the join operation via a clever 
pass through two sorted lists of subset sums. 

In the present paper, we take the meet-in-the-middle approach to counting 
problems, in particular, to counting paths and packings. Here, the join operation 
amounts to consideration of pairs of disjoint subsets of a finite universe, each 
subset weighted by the number of structures that span the subset. We begin in 
Sect. 2 by formalizing this as the Disjoint Sum problem and providing an algo- 
rithm for it based on inclusion exclusion techniques [5-7,17,18,20]. In Sect. 3 
we apply the method to count paths of k edges in a given n-vertex graph in 
time 0*((^'^2))> throughout the paper, O* suppresses a factor polynomial in the 
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mentioned parameters (here, n and k). In Sect. 4 we give another appHcation, to 
count fc-packings in a given family of m-clcment subsets of an n-clcnicnt universe 
in time O* {{jn'k/2)) ■ both problems we also present sHghtly slower algorithms 
that require only polynomial space. 

We note that an earlier report on this work under a different title [8] already 
introduces a somewhat more general technique and an application to counting 
paths. The report has been cited in some recent papers [1, 26], which we, among 
other related previous work, discuss below. 

1.1 Related Work and Discussion 

Deciding whether a given n-vertex graph contains a Hamiltonian path, that is, a 
simple path of n — 1 edges, is well known to be NP-hard. The fastest known algo- 
rithms, due to Bellman [3,4] and, independently. Held and Karp [14], are based 
on dynamic programming across the vertex subsets and run in time 0*(2"). 
Equally fast polynomial-space variants that actually count all Hamiltonian paths 
via inclusion-exclusion were discovered later by Kohn, Gottlieb, and Kohn [20], 
and independently, Karp [17]. Our algorithm (cf. Theorem 2, for A; = n — 1), too, 
runs in time 0*(2"), if allowing exponential space. 

In this light, it is intriguing that the parameterized problem of counting paths 
of k edges seems harder than the corresponding decision problem; this is the 
present understanding that has emerged from a scries of works, starting perhaps 
in Papadimitriou and Yannakakis's [24] conjecture that for k = O(logn) the de- 
cision problem can be solved in polynomial time. The conjecture was proved by 
Alon, Yuster, and Zwick's [2] color-coding technique that gave a randomized al- 
gorithm with expected running time 0*(5.44'^) and a derandomized variant with 
running time 0*(c'^) for a large constant c. With a more efficient color-coding 
scheme, Chen, Lu, Sze, and Zhang [9] improved the latter bound to 0*(12.8'^); 
see also Kneis, Molle, Richter, and Rossmanith [19]. Using completely different 
techniques, Koutis [22], followed by Williams [27], developed a randomized al- 
gorithm that runs in expected time 0*(2'^). Unfortunately, it is unlikely that 
the randomization based techniques extend to counting. For instance, very re- 
cently Alon and Gutner [1] showed that color-coding is doomed to fail as every 
"balanced" family of hash functions from a k-set to an n-set is of size at least 
c(fc)nL'=/2J for some function c. Flum and Grohe [12] proved another negative 
result, namely that the counting problem is #W[l]-hard with respect to the pa- 
rameter k. From a positive side, a very recent result of Vassilevska and Williams 
[26] implies that fc-edge paths can be counted in time O* (2'^(fc/2)!(^"2)) poly- 
nomial space; our polynomial-space algorithm (cf. Theorem 3) is faster still, by 
a factor of (4/3)'=/2(^/2)!. 

Concerning set packings the situation is analogous, albeit the research has 
been somewhat less extensive. Deciding whether a given family of / subsets 
of an n-clement universe contains a fc-packing is known to be W[l]-hard [11], 
and thus it is unlikely that the problem is fixed parameter tractable, that is, 
solvable in time c{k)f'^ for some function c and constant d. If / is fairly large. 
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say exponential in n, the fastest known algorithms actually count the packings 
by employing the inclusion exclusion machinery [5,6] and run in time 0*(2"). 
This bound holds also for the presented algorithm (cf. Theorem 4). 

Again, it is interesting that there is a natural parameterization under which 
counting fc-packing seems harder than the corresponding decision problem. In- 
deed, Jia, Zhang, and Chen [16] showed that the decision problem is fixed pa- 
rameter tractable with respect to the total size mk of the packing, assuming 
each member is of size m. Koutis [21], followed by Chen, Lu, Szc, and Zhang [9], 
gave faster algorithms with running time 0*{c'^'') for some constant c; we note 
that here the running time also grows about linearly in the input size /, which 
can be as large as (^). For counting m-set fc-packings, previous techniques [5, 
6] alone only give a running time bound of if mk < n/2 and 0*(2") 

otherwise. Besides the present work, we arc aware of two recent improvements; 
For the special case of counting t-matchings, that is 2-set t/2-packings,'' Vas- 
silcvska and Williams [26] give a time bound of O* (2*+^^*^ (^"2) ) , where c{t) is 
of the order o{t); our polynomial-space algorithm (cf. Corollary 1) turns out to 
be slightly faster, by a factor of about (4/3)*/^. For the general case, Koutis and 
Williams [23] give a time bound of O* (n™*^/^) ; our bounds (Theorem 4) appear 
to be superior, e.g., when mk grows linearly in n. 

The presented meet-in-the- middle approach resembles the randomized divide- 
and-conqucr technique by Chen, Lu, Szc, and Zhang [9] and the similar divide- 
and-color method by Kneis, Molle, Richter, and Rossmanith [19], designed for 
parameterized decision problems. These can, in turn, be viewed as extensions of 
the recursive partitioning technique of Gurcvich and Shclah [13] for the Hamil- 
tonian Path problem. That said, our contribution is rather in the observation 
that, in the counting context, the join operation can be done efficiently using 
the inclusion-exclusion machinery. While our formalization of the problem as 
the Disjoint Sum problem is new, the solution itself can, in essence, already be 
found in Kennes [18], even though in terms of possibility calculus and without 
the idea of "trimming," that is, restricting the computations to small subsets. 
Kennes 's results were rediscovered in a dual form and extended to accommodate 
trimming in the authors' recent works [5-7]. 

2 The Disjoint Sum Problem 

Given two set famihes A and S, and functions a and (3 that associate each 
member of A and S, respectively, an element from a ring R, the Disjoint Sum 
problem is to find the sum of the products a(A)(3(B) over all disjoint pairs 
of subsets {A, B) in the Cartesian product 7l x S; denote the sum by aM (}. 
In applications, the ring R is typically the set of integers equipped with the 
usual addition and multiplication operation. Note that, had the condition of 
disjointness removed, the problem could be easily solved using about |yi| -I- |!B| 

* Whether counting f-matchings is fixed-parameter tractable remains a major open 
question in parameterized complexity. 
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additions and one multiplication. However, to respect the disjointness condition, 
the straightforward algorithm appears to require about |^||®| ring operations 

and tests of disjointness. 

In many cases, we fortunately can do better by applying the principle of 
inclusion and exclusion. The basic idea is to compute the sum over pairs {A, B) 
with A n B = by subtracting the STim over pairs with A n B = X 7^ from 
the sum over pairs with no constraints. For a precise treatment, it is handy to 
denote by N the union of all the members in the families A and H, and extend 
the functions a and /3 to all subsets of N by letting them evaluate to outside 
A and 23, respectively. We also use the Iverson bracket notation: [P] = 1 if P is 
true, and [P] = otherwise. Now, by elementary manipulation, 

amp = ^ ^ [A n B = 0] a(A) /3(S) 

A B 

A B X 

X A B 

= E(-i)'"'(E"(^))(E/^(^))- (1) 

X ADX BDX 

Here we understand that A, B, and X run through all subsets of N unless other- 
wise specified. Note also that the second equality holds because every nonempty 
set has exactly as many subsets of even size as subsets of odd size. 

To analyze the complexity of evaluating the inclusion-exclusion expression 
(1), we define the lower set of a set family 3^, denoted by J,?', as the family 
consisting of all the sets in 5" and their subsets. We first observe that in (1) it 
suffices to let X run over the intersection of [A and |!B, for any other X has no 
supersets in yi or in B . Second, we observe that the values 

a{X) = E "(^) , 

ADX 

for all X G [A, can be computed in a total of n ring and set operations, as 
follows. Let ai, a2, . . . , a„ be the n elements of N. For any i = 0,1, ... ,n and 
X G lA define ai{X) as the sum of the a{A) over all sets A G lA with ADX 
and An{ai, a2, . . . , Oi} = XD {ai, a2, . . . , Oi}. In particular, q;„(X) = a(X) and 
ao{X) = a{X). Furthermore, by induction on i one can prove the recurrence 

ai-i{X) = [ai ^ X] ai{X) + [XU {ai} e i^l] ai{X U {aj) ; 

for details, see closely related recent work on trimmed zeta transform and Moe- 
bius inversion [6,7]. Thus, for each i, the values ai{X) for all X S J,^ can be 
computed with ||^| ring and set operations. 
We have shown the following. 
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Theorem 1. The Disjoint Sum problem can be solved with 0{n ( \iA\ + )) 
ring and set operations, and with a storage for 0(\iA\ + \i'B\) ring elements, 
where n is the number of distinct elements covered by the members of A and '£>. 

3 Paths 

Consider paths in an undirected graph with vertex set V and edge set E. Define 
a k-edge path as a sequence of A; + 1 distinct vertices vqVi ■ ■ - Vk such that the 
adjacent vertices Vi-i and Vi axe connected by an edge Vi^iVi in E, for i = 
1,2, ... ,k. Wc call the set {vq, Vi, . . . ,Vk} the support of the path and vq and Vk 
the ends of the path. For any vertex and a subset of j vertices S* C y, let Pj{S, v) 
denote the number of j-edge paths with an end v and support S U {v}. Clearly, 
the values can be computed by dynamic programming using the recurrence 

Po{S, v) = [S = 9], pj{S, v) = Y,p,-i{S \ {u}, u) [uv G E] for j > . 
Alternatively, one may use the inclusion-exclusion formula [17, 20] 

res 

where Wj(Y,v) is the number of j-edge walks starting from v and visiting some 
vertices of Y , that is, sequences UqUi ■ ■ • Uj with Uq = v, each Ui^^Ui G E, and 
Ui,U2, . . ■ ,Uj € Y. Note that for any given Y, v, and j, the term Wj{Y,v) can 
be computed in time polynomial in n. Using either of the above two formulas, 
the values Pj{S, v) for all w G V and sets S C V \ {v} of size j, can be computed 
in time here and henceforth, (^^.) denotes the sum of the binomial 

coefficients (g) + (i) H 1" (f) ■ In particular, the number of fc-edge paths in the 

graph is obtained as the sum of Pk{S, v) over all v G V and S CV \ {v} of size 
fc,intimeO*(Q). 

However, meet-in-the-middle yields a much faster algorithm. Assuming for 
simplicity that k is even, the path has a mid- vertex, vi~/2, at which the 
path uniquely decomposes into two fc/2-edge paths, namely vqVi---vi^/2 s-nd 
'fk/2'fk/2+i'''Vk, with almost disjoint supports. Thus, the number of fc-edge 
paths is obtained as the sum of the products 

Pk/2{S,v)pk/2{T,v)/2 

over all vertices v E V and disjoint pairs of subsets S,T C V\ {v} of size k/2. 
Applying Theorem 1, once for each v £ V, with A = 'B = {SCV \ {v} -131 = 
k/2} and a = P = Pk/2 gives the following. 

Theorem 2. The k-edge paths in a given graph on n vertices can be counted in 
time 0*((")). 
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In the remainder of this section we present a polynomial-space variant of the 
above dcsc;ribccl algorithm. Lot the mid-vcrtcx v be fixed. Then the task is to 
compute, for each X <ZV\ {v} of size at most k/2, the sum 

SDK SDXYCS 

in space polynomial in n and k. If done in a straightforward manner, the running 
time, ignoring polynomial factors, becomes proportional to the number of triplets 
{X,S,Y) with X,Y C S C V\{v} and \S\ = k/2. This number is ('^72)2'= 

because there are ('^72) choices for S and for any fixed S, there are 2*^/^ choices 
for X and 2'^/^ choices for Y. 

A faster algorithm is obtained by reversing the order of summation: 



E Pk/2{S,v) = ^«;,/2(y,t;)^(-l)l^\^l [X,YC S] 

SDX Y S 

here Y and S run through all subsets of \ {v} of size at most k/2 and exactly 
k/2, respectively. The latter equality holds bec;ause S is of size k/2 and contains 
X UY. It remains to find in how many ways one can choose the sets X and Y 
such that the union U = X U F is of size at most k/2. This number is 

because there are ("J^) ways to choose U of size s, and one can put each element 
in U either to X or Y or both. 



Theorem 3. The k-edge paths in a given graph on n vertices can be counted in 
time 0*(3'^/^(^"2)) space polynomial in n and k. 



4 Set Packing 

Next, consider packings in a set family J consisting of subsets of a universe N. 
We will assume that each member of 3^ is of size m. A k-packing in is a set of k 
mutually disjoint members of 3^. The members Fi,F2, . . . ,Fk of a fc-packing can 
be ordered in fc! different ways to an ordered k-packing F1F2 • • • Fk- Define the 
support of the ordered fc-packing as the union of its members. For any S C N, 
let TTj{S) denote the number of ordered j-packings in 3^ with support S. The 
values can be computing by dynamic programming using the recurrence 

no{S) = [S = $], 7rj{S)=Y,^j-iiS\F)[Fe3] for j > . 

FCS 
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Alternatively, one may use the inclusion-exclusion formula 

YCS ^ FCY ' 

(here we use the assumption that every member of is of size rri) [5, 6]. Using 
the inclusion-exclusion formula, the values tTj {S) for all 5* C of size mj can be 
computed in time where n is the cardinality of N\ a straightforward 

implementation of the dynamic programming algorithm yields the same bound, 
provided that m is a constant. In particular, the number of fe-packings in J is 
obtained as the sum of TTkiS) / kl over allSCN of size mk, in time O* ((^^J). 

Again, meet-in-the- middle gives a much faster algorithm. Assuming for sim- 
plicity that k is even, we observe that the ordered fc-packing decomposes uniquely 
into two ordered fc/2-packings F1F2 ■ ■ ■ and -F'fc/2+i-Ffc/2+2 " ' Fk with dis- 
joint supports. Thus the number of ordered fc-packings in 3^ is obtained as the 
sum of the products 



7rfc/2(5) 7rfe/2(T)/2 

over all disjoint pairs of subsets S,T C N of size mk/2. Applying Theorem 1 
with A = 'B = {S C N ■.\S\ = mk/2} and a = f3 = nk/2 gives the following. 

Theorem 4. The k-packings in a given family of m-element subsets of an n- 

element set can be counted in time 0*((^^^2))- 

We next present a polynomial-space variant. The task is, in essence, to com- 
pute for each X C N of size at most mk/2 the sum 

E-^m^)=EE(-i)"'"'(E[^^^)''' 

in space polynomial in n, A:, and m. 

As with counting paths in the previous section, a faster than the straightfor- 
ward algorithm is obtained by reversing the order of summation: 



SDX Y ^ FGY ' S 

k/2 



Y ^ FdY ' ^ 



n-\X\JY\ 
mfc/2- |Xuy, 



here Y and S run through all subsets of N of size at most mk/2 and exactly 
mk/2, respectively. It remains to find the number of triplets (X, y, F) satisfying 
|X U y| < mk/2, \F\ = m, and f C F. This number is 



s)\m) 2\mk/2J\ m 

3 / n 
2 \mk/2 



<li T/JS'-'^/'. (2) 
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because there are {^^j^ choices for the union U = Xl^Y oi size s, within which 

there are (j^) choices for F; the elements in F can be put to only Y or to both 
X and Y, whereas each of the remaining s — m elements in U is put to either X 
or Y or both. 

Theorem 5. The k-packings in a given family of m-element subsets of an n- 
element set can be counted in time O* {5"'''/^ {J^^^)) in space polynomial in n, 
k, and m. 

We remark that the upper bound (2) is rather crude for small values of m. 
In particular, provided that m is a constant, we can replace the constant 5 by 3. 

Corollary 1. The k-packings in a given family of 2-element subsets of an n- 
element set can be counted in time 0*(3^(2)) in space polynomial in n and 
k. 
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